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We present a calculation of the correlator (0|T++(r)T++{0)|0> in = 1 SYM 
theory in 2+1 dimensions. In the calculation, we use supersymmetric discrete 
light-cone quantization (SDLCQ), which preserves the supersymmetry at every 
step of the calculation. For small and intermediate r the correlator converges 
rapidly for all couplings. At small r the correlator behaves like as expected. 

At large r the correlator is dominated by the BPS states of the theory. We find a 
critical coupling where the large-r correlator goes to zero; it grows like the square 
root of the transverse cutoff. 



1 SDLCQ of = 1 SYM(2+1) 

In this note we display the results of numerical calculations of the stress-energy correlator 
in A/" = 1 SYM theory in three dimensions. The motivation for this calculation is twofold. 
Firstly, correlators of gauge-invariant operators, such as the stress-energy tensor, are of 
particular interest for string/field theory correspondence scenarios because they can be 
evaluated in both of the (conjectured dual) theories El. Secondly, the correlator uses all 
spectral information, and therefore tests also the wave functions of the states, not only the 
masses. Additionally, it is technically non-trivial to go in a numerical approach to higher 
dimensional problems. _ 

The method of SELLCQ has been described elsewhere 0, so we can be brief here. Following 



where plays the role of a time and x~ is a spatial coordinate. In the sequel, the total lon- 
gitudinal momentum is denoted by , P~ is light-cone energy or the Hamiltonian, and 
the total transverse momentum. It was shown in Ref □ that DLCQ preserves supersymme- 
try, if the supercharge Q~ rather than the Hamiltonian P~ is discretized. Supersymmetry 
is the most important symmetry in this calculation, because it keeps the theory finite even 
at finite discretization parameter, and we shall follow this approach. 

Let us now formulate J\f = 1 SYM(2+1) within this framework. The action is 



where the spinor has two components {ip,x)- We use the light-cone gauge, = 0. We 
express everything in terms of the physical degrees of freedom, namely ip and (p = A-^ , and 
can derive the light-cone supercharges 
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which fulfill the supersymmetry algebra 

{Q+,Q+} = 2\/2P+, {Q-,Q-} = 2V2F-, {Q+, Q^} = -4P_l. 

To evaluate the theory on a computer, we need a finite-dimensional Hamiltonian operator. 
This is achieved by discretizing the theory in the following way. We compactify on a 
circle of period 2L — 2ttK/P^, where the harmonic resolution K is effectively a cutoff in 
particle number. The longitudinal momenta can therefore take integer values rii = 1, 2..., K, 
in proper units. The transverse direction x^ is compactified on circle of period /, with the 
transverse cutoff being 2ttT/1. This cutoff is symmetric, and allows the transverse momenta 
to take values = 0, ±1, ±2, ±T. It is a fundamentally different cutoff in the sense 
that it is only a momentum, and not a particle-number cutoff. 

Close inspection of the resulting supercharges reveals two exact symmetries, namely 
transverse parity, which accounts for a 'parity doubling', and a non-degenerate ^'-Symmetry, 
associated with a flip of the color indices. 



2 Correlation functions 



We want to compute an expression of the form 

F{x+,x-,x-^) = (0|T++(x+,a;",a;^)T++(0,0,0)|0), (1) 
where T++(a;) is a component of the stress-energy tensor, namely the momentum operator 

I 



T++{x) = tr 



^ ^ B 



(x) 



(x) 



Both the boson and the fermion contributions are two-body operators, so that only two- 
particle states will contribute to the correlation function. 

To evaluate the correlator, Eq. (|l|), we take the coUinear limit, x-^ — 0, and insert 
complete set of states |a) with energy eigenvalues P~ into the r.h.s. of Eq. (0). It is 
straightforward to calculate the correlator, Eq. (|l|), in the free case and it then shows the 
expected l/r^ behavior for small distances. The real challenge is the evaluation of the 
correlator when the complete set of bound states is inserted in Eq. (|l|). Wc now have to 
write 



F(x+,x",0) = 



E 

n.m.s.t 



"(n, m)e 



:r++(s,t)|o), (2) 



where Pq„ is Hamiltonian operator. We introduce the convenient notation 



Nu TT 



"(n, to)|0). 



Inserting a complete set of bound states \a) with masses Ma, and evaluating the sums over 
K and N± as integrals yields 



I 



-) Fix\x-,0) 



E 



I 



2(27r)5/2 
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\{u\a)[' 



(3) 



The term in square brackets is constructed in such a way that it is free of any unphysical 
constants. It is calculated numerically, and multiplied by a function including the modified 
Bessel function Kg/2- Collecting powers of r reveals that the individual terms in the sum 
of Eq. (^) behave like in the small r limit. 
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Figure 1. The log-log plot of the correlation function / = {T++ {x)T++ (0)) ^fuT^^i ™- ^ 

for 3=0.10 at K=4; (b) for g=l a,t K = 5. The transverse cutoff runs from T = 1 to 9. 



3 Numerical Results 

Let us now look at the numerical results, Fig. |^. First we focus on the small distance 
behavior of the correlator. Here, the correlator converges from below with increasing T. 
We expect a behavior from the free particle case, but as we saw in the last section, each 
individual bound state behaves like Therefore a coherent behavior of all the states 

has to take place. This makes the reproduction of the 1 /r^ behavior a non-trivial check. It 
turns out that the theory passes this test, see Fig. |l](a): for —0.5 < logr < 0.5 the plot of 
times the correlator falls like l/r. At larger K, Fig. |^(b), we see the same behavior at 
smaller r as expected. 

Now let us look at the correlation fimction at large distances. This region is totally 
determined by massless states. There are actually two types of massless states in the theory. 
The massless states at vanishing coupling g = gyuV^cl/'^'TT^^'^ = are mere reflections of 
the states of the dimensionally reduced theory. They behave as g^Ml_^_l, where Mfj^-^ are 
the masses of the two-dimensional theory □. We therefore expect for g ~ no dependence 
of the correlator on the transverse momentum cutoff T at large r. This is exactly what we 
see in Fig. 0(a) . The other kind of massless states are exactly massless for all couplings and 
are BPS states. Since they are guaranteed to be massless by the BPS symmetry, they have 
to have a complicated dependence on the coupling g through their wavefunction, since we 
see that the large r behavior of the correlator changes with g. In other words, the correlator 
provides us with information on wave functions of the BPS states. 

It turns out that the coupling dependence of the large-r limit of the correlator is very 
interesting. In Fig. |^ we show the large r limit of the correlator as a function of the coupling 
for two different longitudinal cutoffs. Surprisingly, the logarithm of the correlator does not 
change monotonically with but has singularity at a 'critical' coupling which is a function 
of K and T. If we plot the 'critical' couplings vs. ^/T we find that coupling is a linear 
function of ^/T and the coefficient is largely independent of the cutoff K. 
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4 Conclusions 



In this note we showed the first calculation of the stress-energy correlator from first principles 
in three-dimensional Af = 1 SYM theory. We recover the behavior of the free particle 
correlator, because the individual behaviors of the bound states add up coherently. 
The BPS states show an interesting behavior: although their masses are fixed by the BPS 
symmetry, their wavefunctions have to depend on the coupling, as can be inferred by the 
correlator data. Surprisingly, we find a vanishing of the correlator at large distances at a 
critical coupling gcrit oc VT- 

A word on the computer code may be in order. The latest version of the code handles 
two million states, and uses all the present symmetries, which in turn reduces the size of the 
matrix to be diagonalized by a factor of eight. We are therefore optimistic that calculations 
in full 3-1-1 dimensions are within reach. Another interesting theory in three dimensions is 
A/" = (8, 8) SYM(2-|-1), which is conjectured to correspond to a string theory of D2-branes, 
and we are working to generalize our code to tackle the problems associated with a large 
number of 'fiavors'. Our current work should result soon in a complete (low- lying) spectrum 
of A/" = 1 SYM(2-|-1), including all the information in the wavefunctions. These results could 
then be compared to future lattice results, namely the supersymmetric glueball spectrum. 
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